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Lattice QCD with Dynamical Quarks from the UKQCD Collaboration 

Chris Allton 

UKQCD Collaboration 

Department of Physics, University of Wales Swansea, Swansea SA2 8PP, United Kingdom 

A brief overview of the lattice technique of studying QCD is presented. Recent results from the UKQCD 
Collaboration's simulations with dynamical quarks are then presented. In this work, the calculations are all at 
a fixed lattice spacing and volume, but varying sea quark mass from infinite (corresponding to the quenched 
simulation) down to roughly that of the strange quark mass. The main aim of this work is to uncover dynamical 
quark effects from these "matched" ensembles. 



1. Introduction 

In lattice QCD, the quark degrees of freedom 
are defined on the sites of the 4-dimensional lat- 
tice, and the gluonic degrees of freedom are de- 
fined on the links between the lattice vertices. 
This enables a lattice version of the theory to be 
defined which maintains as many of the symme- 
tries of the continuum theory as possible. (For a 
reviews of lattice QCD see jl]].) 

As in other theoretical field theory techniques, 
typical quantities of interest can then be ex- 
pressed as n-point functions, G n (t), of hadronic 
operators. Quantities of physical interest, such as 
hadronic masses, can be extracted from calcula- 
tions of the 2-point function G?2(i) as follows. It 
can be shown that (in Euclidean space-time), 



G a (f) 
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The sum is over all states \i > with the same 
quantum numbers as the operators which define 
the 2-point function, G-2.it)- Hence the ground 
state mass of a hadron, Mq, can be extracted 
simply by fitting G%(t) to an exponential at suf- 
ficiently large times t. 

This is not the end of the story though because 
Mq (and any other quantity extracted from a lat- 
tice calculation) is a function of the input param- 
eters of the simulation. These are the coupling 
strength go, the quark mass m and the volume 
L 3 . Ideally what is required is the continuum 
limit, a — > of Mq (a is the lattice spacing) which 



corresponds to go — > (since QCD is asymptoti- 
cally free). However, as has been found by several 
groups (e.g. tne (dynamical or sea) quark 

mass m sea also has impact on the lattice spacing 
a, i.e. a = f(g ,m sea ). 

Another issue regarding lattice QCD is the er- 
ror which is introduced in the discretisation of the 
continuum theory. In the conventional approach, 
the lattice action reproduces the continuum ac- 
tion only up to terms of 0(a). Over the last 
several years, lattice groups, including UKQCD, 
have been using so-called improved versions of the 
lattice action which have lattice errors of (D(a 2 ). 
While this is an obvious improvement over con- 
ventional lattice actions, it does not preclude lat- 
tice results being contaminated with lattice arte- 
facts, all-be- it at the C(a 2 ) level. 

In order to establish a clean signal it is imper- 
ative to fix a — const for two reasons: (i) since 
the lattice calculation leaves the answer correct 
to 0(a 2 ), to keep this error constant, a must be 
fixed; (ii) in order to maintain fixed finite-size ef- 
fects, a must be fixed (simply because the lattice 
size L = Na where N is the number of lattice 
sites along a side). Since a = f(go,'m sea ) (see 
above) the requirement that a = const defines a 
"matched" trajectory in the (go,m sea ) plane. 

The UKQCD Collaboration has embarked 
upon a study of dynamical quark effects in QCD 
lattice simulations by simulating along a matched 
trajectory. This paper presents an overview of 
this work. See M for a full description. 
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2. Simulation Details 



The standard Wilson gauge action was used 
together with the Clover C(a)-improved fermion 
action. The coefficient, csw, used in the Clover 
term was non-perturbatively determined by the 
Alpha Collaboration [||. 

A range of (f3, K sea ) values were chosen in order 
to maintain a constant value of the lattice spacing 
(/3 = 6/g%). This used the technology of @. The 
parameter values for the simulations are displayed 
in Table |l|. The lattice volume used was 16 3 x 32 
throughout. 

Table [l] lists also the lattice spacing obtained 
from the Sommer scale, tq |Q. The first error 
listed is statistical, and the second (shown as ±*) 
is the systematic error from variations in the fit 
used for the static quark potential. The central 
values quoted were obtained using all potential 
data satisfying \/2 < r < 8. It can be seen 
that the last four simulations listed (i.e. those at 
K soa = 0.1350,0.1345,0.1340 and 0) are matched 
to within errors, with the K sca = simulation 
corresponding to a quenched run. 

The simulations at K soa = 0.1355 and 0.13565 
were performed in order to study lighter sea quark 
masses and do not lie on the above matched tra- 
jectory. 

Full details of autocorrelation times and other 
algorithmic issues can be found in Q . 

3. The static potential 

The static quark potential was determined for 
all our ensembles using the method originally 
proposed in |^] and using a variational basis of 
"fuzzed" links §. 

Figure 1 plots the static potential in units of 
ro . The zero of the potential has been set at r = 
tq. The data are well described by the universal 
bosonic string model potential JTo| which predicts 

[V(r)-V(r )]r = (1.65 - e) (j- - lj 



- 1 



(2) 



Of course, the fact that the scaled potential mea- 
surements all have the same value and slope at 
r = tq simply reflects the definition of tq. Figure 
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Figure 1. The static QCD potential expressed in 
units of tq. The dashed curve is a string model 
described in the text. 



■ 5.93 Quenched 

5.29, c=1.92, k=. 13400 
i 5.26, 0=1.95, k=. 13450 
• 5.20, c=2.02, k- 13500 



Figure 2. The difference between the static QCD 
potential expressed in physical units and the pre- 
diction of the string model described in the text. 
For clarity, only data from the matched ensembles 
are shown. 
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Table 1 

Lattice parameters together with measurements of a and quark masses. 



p 




csw 


#conf. 


a ro [Fm] 


aj [Fm] 


-^unitary j ^unitary 


5.20 


0.13565 


2.0171 


244 


0.0941(8)^ 3 






5.20 


0.1355 


2.0171 


208 


0.0972(8)^ 


0.1101* 


o.578±tl 


5.20 


0.1350 


2.0171 


150 


0.1031(09)1?° 


0.11511 


0.7001;? 


5.26 


0.1345 


1.9497 


101 


0.1041(12)1^ 


0.1181= 


0.7831= 


5.29 


0.1340 


1.9192 


101 


o.ioi8(io)if 


0.11611 


0.8351^ 


5.93 





1.82 


623 


0.1040(03)1* 


0.11861^ 


1 



2 shows the deviations from this model potential 
for the matched ensembles. Here e — tt/12 ||l0| . 

Overall there is no discernable difference be- 
tween the ensembles at distance r ~ r . Fur- 
thermore, the data follow the string model very 
well. However, at shorter distances, r < 0.5ro, 
there is a deviation from the string model, and 
a variation amongst the ensembles. This seems 
to be systematic in the sea quark mass - the de- 
viation from the string model increases as m sea 
decreases. In fact, careful correlated fits of the 
potential show that the parameter e in Eq. (Q) in- 
creases by 18±\l% in going from the quenched to 
K sca _ 23500 data. We note that the ensem- 
bles in Figure 2 are matched, and therefore we 
can exclude lattice systematics from the effects 
that we have observed. (Similar findings in the 
case of two flavours of Wilson fermions have been 
reported by the SESAM-TxL collaboration |I) 
where an increase of 16 — 33% was found.) Fi- 
nally we note that there is no evidence of string 
breaking in our static quark potential. However 
the distance scales covered r 1.3/m is not large. 

4. Hadronic Spectrum 

In this section one of our main aims will be to 
uncover unqucnching effects in the light hadron 
spectrum. Because we have a matched data set, 
any variation amongst our ensembles can be at- 
tributed to unquenching effects. However, the 
task of identifying variations is likely to be hard 
for those quantities which are primarily sensitive 
to physics at the same scale as that used to define 
the matching trajectory in the ((3, re sea ) parame- 
ter space (ro in this case). This is expected to be 



the case for the hadron spectrum considered here 
where the quark masses are still relatively heavy. 

Two-point hadronic correlation functions were 
produced for each of the ensembles using interpo- 
lating operators for pseudoscalar, vector, nucleon 
and delta channels as described in Mesonic 
correlators were constructed using both degener- 
ate and non-degenerate valence quarks, whereas 
only degenerate valence quarks were used for the 
baryonic correlators. Full details of the fitting 
procedure can be found in Q. 

4.1. The J parameter 

In Figures 3 and 4 the vector meson masses and 
hyperfinc splittings are plotted against the corre- 
sponding pseudoscalar masses for all the datasets. 
It is difficult to identify an unquenching signal 
from these plots - the data seem to overlay each 
other. Note that in Q, it was reported that there 
was a tendency for the vector mass to increase 
as the sea quark mass decreases (for fixed pseu- 
doscalar mass). The observations for the present 
matched dataset imply that this may have been 
due to either an O(a) effect (since the dataset in 
j^] was not fully improved at this level) or a finite 
volume effect. The conclusion therefore is that it 
is important to run at a fixed a in order to disen- 
tangle unquenching effects from lattice artifacts 
or finite volume effects. 

A possible explanation for the lack of convinc- 
ing unquenching effects in our meson spectrum 
is the following. Our matched ensembles are de- 
fined to have a common ro value, so any physical 
quantity that is sensitive to this distance scale 
(and the static quark potential itself) will also, 
by definition, be matched. Our mesons, because 
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Figure 3. Vector mass plotted against pseu- 
doscalar mass squared in units of r$, together 
with the experimental data points. 



they are composed of relatively heavy quarks, are 
examples of such quantities, and this is a possi- 
ble reason why there is no significant evidence of 
unquenching effects in the meson spectrum. 

A further point regarding hyperfine splitting in 
Figure 4 is that the lattice data for the matched 
ensembles tends to flatten as the sea quark mass 
decreases. (The quenched data has a distinctly 
negative slope, whereas the ft sca = 0.1350 data 
is flat.) Thus the lattice data is tending to- 
wards the same behaviour as the experimental 
data which lies on a line with positive slope (inde- 
pendent of the value used for ro). This behaviour 
is apparently spoiled by the unmatched run with 
K sea _ 0.1355 (see Figure 4) which has a clear 
negative slope. However, the K sca = 0.1355 data 
does not satisfy the finite volume bound of ||, 
and therefore we can attribute the trend in this 
data to a lattice artefact. 

We now study the J-parameter defined as [O 



J = M 



dM v 



IV 



dM 2 PS 



(3) 



K,K* 



In the context of dynamical fermion simula- 
tions, this parameter can be calculated in two 
ways. The first is to define a partially quenched J 
for each value of the sea quark mass. In this case, 
the derivative in (|3|) is with respect to variations 
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O O 9= 5-20 K = 0.1350 
□ n |i = 5.26 k = 0.1345 
O O f) = 5-29 K = 0.1340 
5.93 Quenched 




Figure 4. The hyperfine splitting. 



in the valence quark mass (with the sea quark 
mass fixed). The second approach is to define J 
along what we will term the 'unitary' trajectory, 
i.e. along K sca = K val . In Table 0, the results 
from both methods are given. These values of 
J are around 25% lower than the experimental 
value J expt = 0.48(2). 



P 


^sca 
fv 


J 


First Approach 


5.2000 


0.1355 


0.32tl 


5.2000 


0.1350 


0.393:™ 


5.2600 


0.1345 


0.3651^ 


5.2900 


0.1340 


0.349t^ 


5.9300 


0.0000 


0.376± 9 12 


Second Approach 






0.351 2 2 



Third Approach 

0.431 2 2 

Table 2 

J values from the various approaches as described 
in the text. 



Finally we note that the physical value of J 
(i.e. that which most closely follows the proce- 
dure used to determine the experimental value of 
Jexpt — 0.48(2)), should be obtained from extrap- 
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olating the results from the first approach to the 
physical sea quark masses. We call this the third 
approach. In order to perform this extrapolation, 
we extrapolate the three matched dynamical J 
values obtained from the first approach linearly 
in (M£f ary ) 2 to (M^ itary ) 2 = 0. Af™ itary is the 
pseudoscalar meson mass at the unitary point. 
The value for J from the third approach is pre- 
sented in Table || and we note that it is approach- 
ing the experimental value for J. 

The results from all three approaches are plot- 
ted in Figure 5, together with the experimental 
result. There is some promising evidence that the 
lattice estimate of J increases towards the experi- 
mental point as the sea quark mass decreases (see 
the J value from approaches 1 and 3). 



Experiment 



Approach 3 ========= 



Approach 2 



O Approach 1 



0.131 0.132 0.133 0.134 0.135 0.136 0.137 0.138 0.139 0.14 



Figure 5. J versus re sea using the approaches 
as described in the text. Note that the quenched 
data points have been plotted at K sca = 0.132 for 
convenience. Approaches 2 & 3 are obtained after 
a chiral extrapolation and are shown as banded 
regions. The experimental value J = 0.48(2) is 
also shown. 



Recently there has been a proposed ansatz for 
the functional form of My as a function of Mp S 
|l4| . Although this ansatz would be interesting 
to pursue, all our data have MpsjMy <; 0.6, and 
for this region, the ansatz of H] is linear to high 
precision. Therefore we have chosen to interpo- 
late our data with a simple linear function and 



await more chiral data before using the ansatz of 
& 

4.2. Lattice spacing 

This subsection presents a determination of 
a from the meson spectrum which complements 
that from tq. 

A common method of determining a from the 
meson spectrum uses the p mass. However, this 
requires the chiral extrapolation of the vector me- 
son mass down to (almost) the chiral limit which, 
as was discussed in the previous subsection, may 
be problematic. An alternative method of ex- 
tracting the lattice spacing using the vector me- 
son mass at the simulated data points (i.e. with- 
out any chiral extrapolation) was given in Esl . 
Using this method, we obtain the lattice spacing 
values as shown in Table [l] labelled a,j . Note that 
these are in general 10-15% larger than the val- 
ues from Section || where the lattice spacing was 
determined from ro- A possible explanation for 
this discrepancy is that the potential and mesonic 
spectrum are contaminated with different 0(a 2 ) 
errors (or that the value tq = 0.49 fm is 10-15% 
too small!). 

In order to investigate unquenching effects in 
the meson spectrum, we define the quantity 



Sij(/3,m s 



1 



a,Q3,m* ea ) 
a,-(/3,ra sea ) ' 



(4) 



where a% is the lattice spacing determined from 
the physical quantity i = {M p , Mk, ■ ■ ■}■ Note 
that when 5ij = 0, the lattice prediction of Mi 
with scale taken from Mj agrees with experiment. 
Thus S is a good parameter to study unquench- 
ing effects. We expect that 8ij((3,m sea ) = 0(a 2 ) 
since we are working with a non-perturbatively 
improved clover action. 



In Figure 6, 



is plotted against 



(aMpj ltaiy )" 2 for the matched datasets. In this 
plot we have fixed j = and the various physical 
quantities i are y/a (the string tension) and the 
hadronic mass pairs (Mk*, Mk) & (M p ,M w ). 
The method used to determine the scale from 
these mass pairs is that of [jl5| . It is worth noting 
that the experimental point on this same plot 
would occur at an x— co-ordinate (depressingly) 



of (aAf™ itary )" 2 = (aMJ 



200. 



l/{aM ps (Li[iitary)]" 

Figure 6. Si as a function of l/(aMp™ tary ) 2 
for i = yfa and the mass pairs (Mk*,Mk) & 
(M p , M„.). Si is defined in eq.(|^) with j = tq. 
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Figure 7. The Edinburgh plot for all the data 
sets. All degenerate K val correlators have been in- 
cluded. The phenomenological curve (from |l7|]) 
has been included as a guide to the eye. 



Figure 6 does not show signs of unquenching 
for quantities involving the hadronic spectrum 
i.e. the mass pairs (M K *,M K ) & (M p ,M„), 
However, there is evidence of unquenching effects 
when comparing the scale from r$ with that from 
y/a since the quenched value of 8^ is distinct 
from the dynamical values. 

4.3. Edinburgh Plot 

In Table |l| the ratios MpsjMy are displayed 



for the unitary case n s 



^val 



As can be seen 



the simulated data is a long way from the experi- 
mental value M^/Mp = 0.18. Figure 7 shows the 
'Edinburgh plot' (Mn I My v.s. M PS /M V ) for all 
the data sets. There is no significant variation 
within the dynamical data as the sea quark mass 
is changed, but the dynamical data does tend to 
lie above the (matched) quenched data. This lat- 
ter feature may be indicative of finite volume ef- 
fects since these are expected to be larger in full 
QCD compared to the quenched case [pJS| . 

5. Chiral Extrapolations 

In j| we used 3 approaches to perform the 
chiral extrapolations of hadron masses: "Pseudo- 
Qucnched" ; "Unitary Trajectory" ; and a "Com- 
bined Chiral Fit" . In this publication we will refer 



only to the last approach. We take 
M(K soa ; K val ) 

= A(n sca ) + B(n sca )M PS (K sca ; n val ) 2 

= A + AiM PS (K sca ; K sca )~ 2 

+ [b + SiAfp S (K soa ; K soa )- 2 ] M PS (K sca ; K va1 ) 2 , 

using the nomenclature M = aM. The first argu- 
ment of M («; sca ; K val ) refers to the sea quark and 
the second to the valence quark. 

The results of these extrapolations are shown 
in Table ||. We stress that this functional form for 
the extrapolation is not motivated by theory, but 
is used as a numerical analysis technique in order 
to test for evidence of unquenching effects. As can 
be seen from Table the parameters A\ and Bi 
are compatible with zero (to 2a) and therefore we 
conclude that there is no evidence of unquenching 
effects (i.e. there is no statistical variation in the 
fit parameters A and B with sea quark mass K sea ). 

6. Conclusions 

This paper attempts to uncover unquenching 
effects in the dynamical lattice QCD simulations 
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hadron 



A 



B 



Vector 

Nucleoli 

Delta 



.4921^ 
.663±ii 
.84± 2 2 



-0.00411 0.61+2 .0151? 
0.006^ 1.23±| -.00111 



-0.002j 



0.91" 1 



•021 2 2 



Table 3 

Fit parameters from the Chiral Extrapolations 



at a fixed (matched) lattice spacing (and vol- 
ume) and various dynamical quark masses. This 
approach allows a more controlled study of un- 
qucnching effects without the possible entangle- 
ment of lattice and unquenching systematics. 

We see some sign of unquenching effects in 
the static quark potential at short distance, but 
no significant sign of unquenching effects in the 
meson spectrum. This is presumably since our 
dynamical quarks are relatively massive, and so 
the meson spectrum is dominated by the static 
quark potential. This potential is, by definition, 
matched amongst our ensembles at the hadronic 
length scale ro, and so any variation of the me- 
son spectrum within our matched ensemble must 
surely be a "higher" order unquenching effect 
which is beyond our present statistics. 

It is interesting to note that other work (using 
the same ensembles) has shown interesting un- 
quenching effects in the glueball and topological 
sector 0]. 
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